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Abstract 

For the heat equation in a bounded domain we give a stability result for a smooth diffusion coeffi- 
cient. The key ingredients are a global Carleman-type estimate, a Poincare-type estimate and an energy 
estimate with a single observation acting on a part of the boundary. 



1 Introduction 

This paper is devoted to the identification of the diffusion coefficient in the heat equation using the least 
number of observations as possible. 

Let fl C M" be a bounded domain of M" with n < 3, (the assumption n < 3 is necessary in order to obtain 
the appropriate regularity for the solution using classical Sobolev embedding, see Brezis [3]). We denote 
r = dfl assumed to be of class C^. We denote by i' the outward unit normal to on T = d^l. Let T > 
and to e (0, T). We shall use the following notations Qo = Q x {0,T), Q = Q x {to,T), Y. ^ T x {to,T) and 
So = r X (0,T). We consider the following heat equation: 

{dtq = V ■ (c(x)V(?) in Qo, 
q{t,x) = g{t,x) on Sq, 
q{0, x) = qn in Q. 

Our problem can be stated as follows: 
Inverse Problem 

Is it possible to determine the coefficient c{x) from the following measurements: 

dAdtq)\ito,T).ro and V(Ag(T',.)),Ag(T',-), 9(T', •) in for T' = 
where Fq is a part of the boundary F of ? 

Let q (resp. q) be solution of p.ip associated to (c, g, qo) (resp. (c,g, qo)), we assume 

Assumption 1.1. • qo belongs to H^{Vl))and g is sufficiently regular (e.g. 3 e > such that g E 

H\o,T,H^/^+%dn))nH^{o,T,H'^/^+%dn))) 
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• c, ceC^in), 

• There exist a constant r > 0, such that qo > r and g > r. 

Note that the first item of the previous assumptions implies that (jl.ip admits a solution in H^{tQ,T, 
(see Lions [E])- We will later use this regularity result. The two last items allows us to state that the 
function u satisfies \Aq{x,T')\ > r > and \Wq{x,T')\ > r > in fl (see Pazy [TS], Benabdallah, Gaitan 
and Le Rousseau [1]). 

We assume that we can measure both the normal flux du{dtq) on Fq C dVl in the time interval {to,T) for 

some to e (0,T) and V(Ag), Aq and Vq at time T' G ito,T). 

Our main result is a stability result for the coefhcient c(x): 

For go in H^{^) there exists a constant C — C{il, F, to, T,r) > such that 

< C\9Adtq) ~ d^{dtq)\l2(^(^t^^T)xro) 
+ C|V(Aq(r',-))-V(Ag(r',-))li.(o) 

+ C\Aq{T', .) - Aq{T', + C\Wq{T', ■) - Vq{T', Oli^jo). 

The key ingredients to this stability result are a global Carleman-type estimate, a Poincare-type estimate 
and an energy estimate. We use the classical Carleman estimate with one observation on the boundary for 
the heat equation obtained in Fernandez-Cara and Guerrero [5], Fursikov and Imanuvilov [S]. Following the 
method developed by Imanuvilov, Isakov and Yamamoto for the Lame system in Imanuvilov, Isakov and 
Yamamoto [11] . we give a Poincare-type estimate. Then, we prove an energy estimate. Such energy estimate 
has been proved in Lasiecka, Triggiani ang Zhang [13] for the Schrodinger operator in a bounded domain in 
order to obtain a controllability result and in Cristofol, Cardoulis and Gaitan for the Schrodinger operator 
in a unbounded domain in order to obtain a stability result. Then using these estimates, we give a stability 
and uniqueness result for the diffusion coefhcient c{x). In the perspective of numerical reconstruction, such 
problems are ill-posed and stability results are thus of importance. 

In the stationnary case, the inverse conductivity problem has been studied by several authors. There are 
different approaches. For the two dimensional case, Nachman |14| proved an uniqueness result for the 
diffusion coefficient c G C^(r2) and Astala and Paivarinta |T] for c G L°°{n) with many measurements from 
the whole boundary. In the three dimensional case, with the use of complex exponentially solutions, Faddeev 
[7], Calderon [S], Sylvester and Uhlmann [TB] showed uniqueness for the diffusion coefficient. 
There are few results on Lipschitz stability for parabolic equations, we can cite Imanuvilov and Yamamoto 
|10j . Benabdallah, Gaitan and Le Rousseau [4]. In [4], the authors prove a Lipschitz stability result for the 
determination of a piecewise-constant diffusion coefRcient. For smooth coefficients in the principal part of 
a parabolic equation. Yuan and Yamamoto |17| give a Lipschitz stability result with multiple observations. 
This paper is an improvement of the simple case in |17j where we consider that the diffusion coefficient is a 
real valued function and not a n x n- matrix. Indeed, in this case, with the method developped by |17j . they 
need two observations in order to obtain an estimation of the i/^-norm of the diffusion coefRcient. In this 
case, we need only one observation. 

Our paper is organized as follows. In Section 2, we recall the global Carleman estimate for p.l[) with 
one observation on the boundary. Then we prove a Poincare-type estimate for the coefficient c{x) and an 
energy estimate. In Section 3, using the previous results, we establish a stability estimate for the coefficient 
c{x) when one of the solutions q is in a particular class of solutions with some regularity and "positivity" 
properties. 
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2 Some Usefull Estimates 



2.1 Global Carleman Estimate 

Wc recall here a Carleman-type estimate with a single observation acting on a part Tq of the boundary F of 
il in the right-hand side of the estimate (see [5]), Let us introduce the following notations: 

let /3 be a C'^{fl) positive function such that there exists a positive constant Co which satisfies 
Assumption 2.1. |V/3| > Co > in H, 9^/3 < on T \ Tq, 

Then, we define (3 ^ (3 + K with if = m||/3||oo and m > 1. For A > and t e {to,T), we define the weight 
functions 

gA/3(x) g2AK _ gA/3(x) 



{t-to){T-t)' " ' (t-io)(r-i)' 
If we set — e^^''g, we also introduce the following operators 

MiiP = \7 ■{cV^)+s^X'^c\VP\^ip^ip + s{dtr])4>, 
M21P = StV' - +2sXipc\7p.V'ip - 2sAVc|V/3pV'- 

Then the following result holds (see [S], 0) 

Theorem 2.2. There exist Xq ~ Xo{Cl,To) > 1, so = •so(Ao,T) > 1 and a positive constant C — C{^1,Tq,T) 
such that, for any A > Aq and any s > sq, the following inequality holds: 

(2.2) l|A^i(e-^''g)|li.(Q) + \\M2{e-^\)\\h^Qj 

+sX^ [[ e^^'^VlVgp dx dt + s^X^ [[ e^^^V^i^p 



< C 



■ T 

sX 



f f e-^'^^ifild^qf dxdt+ ff e-2«'? \dtq~V ■ {cVq)f dx dt 
JtoJro JJq 



for allq e H^{to,T,H^{n)) with q ^ onH. 

2.2 Poincare-type estimate 

We consider the solutions q and q to the following systems 

{dtq^V ■ {c{x)Vq) in Qo, 

q{t,x) ^ g{t,x) on Eo, 

q{0, x) = go in il, 

and 

5tg = V • (c(x)Vq) in Qo, 

(2.4) { qit.x) ^ g{t,x) on Eo, 

(7(0, x) — qo in Q. 



3 



We set u — q — q, y — dtu and 7 = c — c. Then y is solution to the following problem 

r dty = V ■{c{x)Vy)+V ■{^{x)V{dtq)) in Qo, 

(2.5) <^ y{t,x) = on So, 

[ 2/(0,a;) = V-(7(.t)V(<jo(x))), in n. 

Note that with ()2.3p and (|2.4p we can determine y{T',x) and we obtain 

(2.6) y(r', x) = V ■ (7(.t) V(g(T', x))) + V • (c(a;)V(u(T', x))). 

We use a lemma proved in [llj for Lame system in bounded domains: 
Lemma 2.3. We consider the first order partial differential operator 

Pog := Vqo ■ Vg 

where go satisfies 
|V/3-Vgo| ^0. 

Then there exists positive constant, si > and C — C{X,T') such that for all s > si 

^2y2 f e-2s^(T')^(j,,)|^|2 dx dy < C ( g-^.^*^') ^ (T') IP0.9I' dx dy 
Jn Jn 

with T' = T]{T') := ri{x,T'), ip{T') ipix,T') and for g G H^{n). 

We assume 

Assumption 2.4. |V/3 ■ Vq(T')| ^ 0, 

Proposition 2.5. Let q be solution of (2^. We assume that Assumption \2.4\ are satisfied. Then there 
exists a positive constant C — C(T', A) such that for s large enough ( s > si), the following estimate hold 
true 

^2y2 f e-2s^(i^')^(r')(|v7|2 + \j\^)dx<C [ e-2-''(^'V"'(T')(|Vj/(T')r + |2;(T')|') dx 
Jn Jn 

+C f e~2-KT')(|v(AM(r')|' + |AM(T')|' + VM(Tf) dx 
Jn 

for-i&HliSl). 

Proof. We are dealing with the following first order partial differential operators given by the equation (|2.6p 

n 

^^0(7) := E d.MT')d..l = y{T') - jAq{T') - \/{cVu){T'). 
We apply the lemma 12.31 for this operator and we can write : 
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s2A2/^e-W')<^(r/)|^|2 < (^J^g-2s„(T')^-l(J^/) (|y(I./)|2_^ |^|2) 

(2.7) 

+C^e-2.^(^') (|Au(r')l' + |Vu(r')l') f^a; 

In the other hand, we use the Xj-derivative of the previous equation (|2.6p . So, for each j we deal with the 
foUowing first order partial differential operator : 

Then under assumption (|2.4p : 



s 



n 



+c/ e-2^''(^')^-i(r')(|a.,7p + l7p + |V7p + |a.,^^p) dx 



n 



So, adding for all j, we can write 



(2.8) s^X^ / e-2.^(^')(/p(r')|V7p dx < C e-^'^^^"^^-\T')\VyiT')\^ dx 

Jn Jn 

+C f e-2-'(T')^-i(2^')(|v^|2 + |^|2 + |V(AM(T')|' + |Aii(T')n dx 



Taking into account (|2.7p and (|2.8p and for s large enough, we can conclude. □ 
2.3 Estimation of f^e-^'''^'^'^\y{T')\'^ dx 

Let r' = i(r + to) the point for which $(t) = has its minimum value. 

We set tp = e~'*''2/. With the operator 

(2.9) M2V' = dti' - +2sXipcVp.Vip - 2sAVc|V/3p?/', 

we introduce, following [2], 



rT' 

1^1 I A'h'ipip dxdt 



to Jo, 



We have the following estimates. 
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Lemma 2.6. Let A > Ai, s > si and let a, b, c, d G L°°{fl). Furthermore, we assume that uq, vq in 
and the assumption is satisfied. Then there exists a constant C ~ C{VI^IjJ^T) such that 



(2.10) 



-2s„(T',.) \y{T',x)\'' dX<C 



to Jt, 



-2s?) 



(p\d^y\^ dx dt 



T 

to Jn 



+S-1/2A-1/2 / / ^-2sn (|^|2 ^ |y^|2) 



Proof. If we compute X, we obtain : 



-2sr]{T',x) 



T' 



yiT',x)f dx = -21 

■T' 



-4sA / / ip cVP -V^jj ^ dx dt- 4sA^ / / ip c|V/3nV'| dx dt. 
Jto Jfi Jto Jn 

Then with the Caiicman estimate ()2.2p . we can estimate all the terms in the right hand side of the previous 

equality and wc have 

-2sr,(T',.) \y^T\x)\Ux<Cs-^/^X-'- l^\Ahn^ + s^X' jj e^'^V'lyP dx dt)j 
+Cs-^\~^'^ (s\ f f e-2^'V |Vyp dx dt + s^X^ [ [ e-^^V^lyP dx dt 



Q 



+Cs-^\-^ /a4 



^''V^lyP dx dt 



Finally, we obtain 



n JtoJTo JJq 

where f — V ■ {jVdtq). Wc assume that g is sufficiently smooth in order to have Vdtq and Adtq in 
L^{0,T,L°°{n)). 

Moreover taking into account that e"^'"''-*-' < e^^*''^^ \ the proof of Lemma is complete. □ 



2.4 Estimation of e-2^^(^'V"Hr')|V?/(T')P 



We introduce 
(2.11) 

In this section, we give an estimation for the energy E{t) at T' . 



E{t)= / c (/j~i(a;,t)e-2-''(^'*)|Vy(x,t)p dx. 



Theorem 2.7. We assume that Assumptions [7771 are checked, then there exist Ai = Ai(f2,a;) > 1, si = 
si(Ai,r) > 1 and a positive constant C = C(r2, Fq, Co, T) such that, for any A > Ai and any s > si, the 
following inequality holds: 



(2.12) E{T') < C 



sA 



[ [ e-2"'V|9^yP dx dt + s [[ 
Jto Jto JJq 



>-2s»?n,,|2 



71 + Nil ) dx dt 
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Proof. We note / = V • {-/{x)\/dtq). 

We multiply the first equation of (j2.5p by e~^'"'V • {c'Vy)Lp~^ and integrate over (ioj?^) x we have : 

(2.13) [ [ ip-^e-^'^^S/ ■ {c\7y)dty dx dt = [ [ ^-^e-^'^^lV ■ (cVy)p dx dt 
J to Jn J to Jn 

e-^'^^^-^V ■ {cVy)f dx dt. 

ito Jn 

we denote A :^ ( ( e'^^'V^^V • {cVy)dty dx dt. 
J to Jn 

Integrating by parts A with respect to the space variable, we obtain 

(2.14) A= [ [ ce-^'^^^-^VydtiVy) dx dt + 2sX [ [ c e-^'^^S/ydtyS/p dx dt 

J to J ft J to Jfi 

-X [ [ c e~2«')(p-iVyat?/V^ dx dt. 

J t o J ft 

Observe that 



Hence, the first integral of the right-hand side of (|2.14p can be written as 

'' c e-2''V"^Vyat(V?/) dxdt^ I [ c e-''''ip~i\/ydt{Vy)e-''''ip-i dx dt 



= / [ ce~'"^>f-i\7ydt{e~'''^ip-^'^y) dx dt + s [ j c e~'^''^^-^\Wy\^dtTi dx dt 
J to Jn J to Jfi 

(2.15) + - / / c e-^'^^ip-^lVyl^dtip dx dt. 

2 J to Jn 

Using an integration by parts with respect the time variable, the first term of (j2.15p is exactly equal 
^E{T'), since E{to) = 0. Therfore, the equations H^J^ . ([^1^ and (|TTS)) yield 

E{T') = -2s [ [ c e~2"'V~^|V2/pc't?7 dx dt - [ [ c e^^'^^ip-^lVyl^dtip dx dt 
J to Jn J to Jn 

-4sA j j c e-^'^^'VydtyVP dxdt + 2\ [ [ c e-^'^^ip-^VydtyVP dx dt 
J to Jn Jto Jn 



+2 f f ip~^e~^''^\V ■ {cVy)\^ dx dt + 2 f f e-^'^^cp-^V ■ {cVy)f dx dt 
Jto Jn Jto Jn 

(2.16) =/l+/2+/3 + /4 + /5+/6- 
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Now, in order to obtain an estimation to E{T'), we must estimate all the integrals It, 1 < i < 6. 

Using the fact that \dtri\ < C{il,uj)T(p'^ , we obtain, in first step, for the integral /i, the following estimation 

\h\<Cs j [ c e^^^'VlVyp dx dt 
J to Jn 



sX^ I I e-^'^^^lVyp dx dt 
Q 



In a second step, the Carleman estimate yields 



\h\<cx- 



sX 



[ j e-^''^ip\d^y\^ dx dt+ f f e~^''^\f\^\ dx dt 

JtaJTo JJq 



where C is a generic constant which depends on fi, Fq, Cmax and T . 
As the same way, we have, for /2, the following estimate 



I/2I <Cs-^X- 



sX / e-^''^ip\d^y\'^ dx dt+ // e-'^'''^\f\'^\ dx dt 



The last inequality holds throught the Carleman estimate and the following inequality 



\dM<c{n,To)T 



.3 V 



Using Young inequality, we estimate /s. 
We have 



\h\ < Cs 



e'^'^l/n dx dt 



< Cs 



sX I I e-'^'^ifild^yl'^ dx dt 



For the integral /4, we have 

I/4I <c 



sX^ II e-^'''^ip-^\yy\'^ dx dt + s-^ 
Q 

T 

<C sA 



-2sri -1 



if \dty\ dx dt 



[ [ e-^'Wd^l^ dx dt + [ [ e-^^^l/H 
JtrJr„ J Jo 



where we have used, for the term containing |Vyp, the following estimate 



16 



we have immcdiatly the following estimate for I5 



1/5 1 < Cs 



sX 



-2S7] 



p\d^y\^ dxdt+ // e-2""|/n dx dt 
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Finally, for the last term /g, we have 



Ihl < c 



|V • (cV2/)p dxdt + s / / e^^-^^l/l^ dx dt 



< C 



sX I I e-^'^^^ld^yl^ dx dt + s 1 1 e-'^''''\f\'^\ dx dt 



The last inequality holds using the following estimate 

ip-^ < c{n,uj)T' 



2V 



If we come back to ()2.16p . using the estimations of J^, 1 < i < 6 and expanding the term /, this conlude the 
proof of Theorem 12.71 □ 



3 Stability Result 

Theorem 3.1. Let q and q be solutions of i2.3\) and {2-4^ such that c — c E Hq{VI). We assume that 
Assumptions [TTTl are satisfied. Then there exists a positive constant C = C(r2,ro,T) such that for s and A 
large enough, 



'l\ p-2s7)(T') 



|c-2^2^|V(c-2)r) 



dx dy<C f [if e-^'^^d^P \d^{dtq - dtq)\^ da dt 



+C / e-2-''(^')(|V(Au(T')P + |Au(T')|' + Vu(T'|2) dx 



Proof. Using the estimates (|2.12p . (|2.10p and Proposition (|2.5p . we obtain 

^2^2 f e-2-''(^')(^(T')(|V7|2 +|7|2) dx<C [ e'^' V~'(r') (|Vy(r')P + ly(r')l') dx 



< C 



+C f e-^-''^^') (|V(Au(T')|' + |Au(r')l' + Vu(T'|2) dx 
Jn 

sX f f e~2"'V|9i.2/P dxdt + s ff e-^'^^djl^ + | V7P) dx dt 



/ / e-^^VI^.yP dx dt + s-^I^X-^'^ / / e-^--? (|^|2 + |V7|2) dx dt 
Jto JTo Jto Jn 



+C [ e-^'"'^^') (|V(Aw(r')P + |Au(T')P + Vu(r'p) dx. 
Jn 

So we get for s sufficiently large 

s^X^ f e-2^''(^'V(r')(|V7p + |7n dx < CsX C f e'^''' ^\d,y\^ dx dt 
Jn Jto Jto 



'to JTo 

+C [ e-2'"'(^') (|V(Aw(r')P + |Au(T')P + Vu(r'p) dx, 
Jn 



and the the theorem is proved. 



□ 



Remark 



• All the previous results are available for fl C M" be a bounded domain of M" with n > 3 if we adapt 
the regularity properties of the initial and boundary data. 

• We give a stability result for two linked coefficient (c and Vc) with one observation. Note that for two 
independent coefficients, there is no result in the litterature with only one observation. 
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